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ABSTRACT 

We report the results of a measurement of susceptibilities in noncompact 
QEDi in 8^ 10"^ and 12^ lattices. Due to the potentiahties of the MFA 
approach, we have done simulations in the chiral limit which are therefore free 
from arbitrary mass extrapolations. Our results in the Coulomb phase show 
unambiguously that the susceptibility critical exponent 7 = 1 independently 
of the flavour symmetry group. The critical couplings extracted from these 
calculations are in perfect agreement with previous determinations based on 
the fermion effective action and plaquette energy, and outside the predictions 
of a logarithmically improved scalar mean field theory by eight standard 
deviations. 
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Non compact Electrodynamics in four dimensions {QED4) has a strongly 
coupled phase where chiral symmetry is spontaneously broken. The strong 
coupling phase is connected to the weak coupling phase through a continuous 
phase transition characterized by a critical point where a continuum limit can 
be defined. The possible existence of a non trivial continuum limit for this 
model is an appealing question which has attracted much attention in the 
last years and whose answer is not well known yet. 

The main difficulty in establishing the nature of the continuum limit of 
this model comes from the fact that it can not be investigated in pertur- 
bation theory. The possible existence of a non gaussian fixed point implies 
that operators of dimension higher than four could become nonperturba- 
tively renormalizable and therefore dimensional analysis does not help us 
to establish what the relevant couphngs are. 

Due to the high complexity of this problem, the first approaches have been 
concentrated in the determination of thermodynamical quantities, critical 
exponents and their comparison with the predictions of a mean field theory 
plus logarithmic corrections. Several numerical methods have been used by 
different groups to simulate dynamical fermions. In refs. [0, the hybrid 
Monte Carlo algorithm was used whereas in our previous work the MFA 
approach ||5| was employed. Even if the numerical results obtained by all 
these groups are in reasonable agreement, their physical interpretations are 
in conflict. 

In Pi the value of the critical coupling is left as a free parameter in a flt 
of the numerical results with a logarithmically improved scalar mean fleld 
theory. The flve parameters flt accommodates quite well the experimental 
data and the prediction of the flt for the critical coupling Pc in the four flavour 
theory is = 0.186(1). These results have been interpreted by the authors 
of 1^ as strong evidence of a trivial continuum limit for this model. Notice 
however that if as recently argued 0, triviality manifests in a different way 
in scalar and fermionic models, the logarithmically improved mean fleld test 
of would not be justifled. 

Conversely, the results reported in and |^ for the critical couplings, 
both of them in very good agreement, show a value for (3c larger than 0.186. 
The excellent agreement between the results of 0] and should be of great 
physical relevance since both simulations are perfectly uncorrelated. In fact 
these two groups use different approaches to simulate dynamical fermions, 
different lattice sizes, different operators and different bare fermion masses 
(the simulations in were done in the chiral limit). 

Since the numerical determination of critical exponents is very sensitive to 
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the value of the critical coupling Pc and that values of (3c > 0.186 could pro- 
duce critical exponents outside of the mean field range, it is very important 
to have more independent measurements of f3c- 

We report here the results of a measurement of the longitudinal, trans- 
verse and non linear susceptibilities in QED4 with 0, 2 and 4 flavours. This 
is the first time that susceptibilities in noncompact QED are calculated in 
the chiral limit (Coulomb phase) and the results are certainly encouraging. 
In fact the values obtained for the critical couplings are in perfect agree- 
ment with previous determinations based on measurements of the fermionic 
effective action and plaquette energy |Q and therefore we confirm again from 
these results that (3c in the four flavour model is definitely larger than 0.186. 
These measurements show also unambiguously that the susceptibility critical 
exponent 7 is 1, independently of the flavour symmetry group. Furthermore 
the value of (3c in the quenched model extracted from the susceptibility is also 
in perfect agreement with previous determinations of it from the computation 
of the chiral condensate 0] and larger than the value (3c = 0.244 in the zero 
flavour model extracted from the singular behavior of the fermionic action 
[|], the last being also in perfect agreement with the critical (3 obtained from 
the monopole percolation transition 0. 

We have done calculations in 8^, 10^ and 12^ lattices and dynamical 
fermions are simulated by mean of the MFA method, the only available nu- 
merical approach which allows to perform realistic numerical simulations in 
the chiral limit. MFA (Microcanonical Fermionic Average) has been exten- 
sively discussed and tested in the literature. Let us however remember 
what are the main ingredients of the MFA approach. 

The basic idea in MFA is the exploitation of the physical equivalence 
between the canonical and the microcanonical formalism, in our case the 
introduction of an explicit dependence on the energy in the computation of 
the partition function. Indeed it can be written as follows: 

Z{(3,m) = J dEn{E)e^^'^^MA{E,m) (1) 

where 

n{E) = J [dA^]Si6VE - Sg[A,]) (2) 
is the density of states at fixed energy E and 
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is the fermionic determinant averaged over field configurations of fixed energy 
E. Note that detA(£',m) does not depend on p. The fermionic determinant 
for a single configuration is a polynomial in the mass: 

detA[A^,m] = J2Cn[A^]m'' (4) 

n 

A modified Lanczos algorithm [^] is used in order to obtain the complete 
set of eigenvalues of the massless fermion matrix, and then we can reconstruct 
the fermionic determinant at any value of the fermion mass m. 

The most general expression for the vacuum expectation value of any 
operator O, after integration of the Grassmann variables, can be written as 

/ n >= ^ ' detA lr\ 

where S^^^iE^m) in (5) is the fermion effective action defined as 

S^.AE, m) = - log d^A(^, m) (6) 



and OdetA means the mean value of the product of the operator O times 
the fermionic determinant, computed over gauge field configurations at fixed 
energy E. 

Since we are interested here in the computation of susceptibilities in the 
chiral limit, we will use the previous expression for the particular cases in 
which O is the longitudinal, transverse and non linear susceptibility. 

In the Coulomb phase, characterized by a ground state invariant under 
chiral transformations in the chiral limit, the longitudinal and transverse 
susceptibilities are equal except a sign. They can be computed by taking for 
the operator O the expression 

where the sum in (7) runs over all positive eigenvalues of the massless Dirac 
operator. 

The transverse susceptibility xt in the broken phase diverges always be- 
cause of the Goldstone boson. The longitudinal susceptibility xl on the other 
hand, can not be computed in the broken phase using equation (7) since in 
this phase the ground state is not invariant under chiral transformations. In 
fact if we work directly in the chiral limit, we average over all ground states. 
In other words there are extra contributions to Xl in the massive case, which 
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vanish in the chiral hmit in the Coulomb phase because of the reahzation of 
chiral symmetry, but giving an important contribution to Xl in the broken 
phase. The standard procedure in the broken phase is then to compute the 
susceptibihty at non zero fermion mass and to extract its value in the chiral 
limit by fermion mass extrapolations. 

Since the susceptibilities in the symmetric phase are free from mass ex- 
trapolations, we will devote the largest part of this paper to the analysis of 
our results in the Coulomb phase which are the more relevant from a physical 
point of view. 

The results in the 8^ lattice are the best from a statistical point of view 
since in this lattice we have diagonalized 3000 completely uncorrelated gauge 
configurations for each value of the energy and repeated this process for 17 
different energies between 0.5 and 1.2. In the lO'' lattice, 300 — 400 gauge 
configurations were diagonalized at each energy, whereas our results in the 
12^ lattice are statistically the poorest. 

In Fig. 1 we have plotted the inverse susceptibility in the chiral limit 
against 13 for the four-flavour model in a 10^ lattice. The continuous line is 
a fit of all the points for 13 > 0.210 with a function 

v-i = ^ ~ ^) rs^ 

^ + r) (/? + r) 

The use of this fitting function is strongly suggested by the results of 
the inverse susceptibility as a function of the energy (Fig. 2), which are 
extremely well fitted in the Coulomb phase by a straight line, and by the 
relation between the mean plaquette energy and P in the massless case 
which is extremely well reproduced by the function 

the value of r depending on the flavour number. 

Fig. 2 contains the results for the inverse susceptibility against the 
energy E in the four and zero flavour models. Again in these cases, all the 
points in the Coulomb phase are very well fitted by a straight line, which 
implies 7 = 1 independently of the flavour number. 

The results of these fits for the critical coupling in the four flavours model 
are Pc = 0.202(3)(8^), pc = 0.2026(21)(10^). The ^ are 2.8, 0.56 respec- 
tively, which show the high reliability of these fits as well as of the value 
7 = 1 for the susceptibihty critical exponent. 
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The previous values of /3c in the 8^ and 10^ lattices are on the other hand 
in perfect agreement with the corresponding values extracted, following the 
method reported in |Q, from the singular behavior of the fermionic effective 
action (see also Table I). As anticipated in the introduction, these values are 
out of the value 0.186(1) reported in by eight standard deviations. 

The values of the critical couplings and energies for different lattice sizes 
and flavour number are reported in Table I, where we include also the critical 
values obtained from the singular behavior of the fermion effective action. As 
can be seen in this Table, the results for the critical couplings extracted from 
the susceptibility and fermion effective action in the two flavour model are 
also in very good agreement, but they are incompatible in the quenched case. 
Notice however that the quenched value of /3c obtained from the susceptibility 
is the same as the one obtained in an independent calculation of the chiral 
order parameter in 24^ lattices [0] whereas the value obtained from the 
anomalous behavior of the fermionic action agrees very well with the critical 
coupling for the monopole percolation transition in the quenched model 



It is not clear to us the physical origin of the discrepancy between the crit- 
ical couplings obtained by measuring different operators in the zero flavour 
model. However the agreement between critical couplings for monopole per- 
colation and singular behavior of the fermionic action is, we believe, a clear 
signal of strong correlation between these two phenomena. 

We have also done measurements of the longitudinal susceptibility xl 
for massive fermions as well as of the non linear susceptibility Xni, which is 
related to the four point function and therefore to the low energy renormal- 
ized coupling. A detailed analysis of these calculations will be reported in 
a longer publication. However we report in Fig. 3 some results obtained 
for Xl^ as a function of the gauge coupling P at two values of the fermion 
mass m = 0.00625,0.01. The value of the dynamical fermion mass in this 
figure (the mass in the integration measure) has been taken equal to zero. 
The inverse chiral susceptibility shows a minimum pointing to the critical /3 
previously obtained from the fit of the results for the massless susceptibility 
in the Coulomb phase. 

Concerning the non linear susceptibility Xni, we would like to anticipate 
that the results for Xni suffer, as expected, from strong finite size effects [p]] . 
One of the contributions to the massless non linear susceptibility operator 
can be defined as 

o„,4e^ (10) 



6 



where again in this case the sum runs over all positive eigenvalues of the 
massless fermion matrix. 

The expression given by equation (10) diverges in the thermodynamical 
limit due to the presence of zero modes both in the Coulomb and broken 
phases. However the divergency in the broken phase is much faster because 
chiral symmetry is spontaneously broken in this phase. Then a rough esti- 
mate of the critical coupling can be obtained from the results for the v.e.v. 
of the inverse of this operator, by defining it as the point where the slope 
of this v.e.v. as a function of the gauge coupling /? (see Fig. 4) takes its 
maximum value. This criterion gives for Pc in the quenched model a value 
between 0.25 — 0.26 in good agreement with the results reported in Table 1. 

In conclusion we have shown how the analysis of the massless chiral sus- 
ceptibility in the Coulomb phase of noncompact QED, even in rather small 
lattices as ours, can be very useful to get precise determinations of the criti- 
cal couplings. Indeed what gives high reliability to the susceptibility fits we 
have done in the Coulomb phase is the fact that all the points, except those 
very near to (3c which could be affected by finite size effects, are very well 
fitted by a single power law term without higher power corrections, even far 
from /3c. 

The numerical simulations quoted above have been done using the Trans- 
puter Networks of the Theoretical Group of the Frascati National Labora- 
tories and the Reconfigurable Transputer Network (RTN), a 64 Transputers 
array, of the University of Zaragoza. 

This work has been partly supported through a CICYT (Spain) - INFN 
(Italy) collaboration. 
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Figure captions 



Fig. 1. Inverse susceptibility in the chiral limit against (3 for the four flavour 
model, 10*^ lattice. 

Fig. 2. Inverse susceptibility as a function of the energy E 'm. & 10^ lattice 
for zero and four dynamical flavours. 

Fig. 3. Inverse longitudinal susceptibility for m — 0.01 and m — 0.00625 

against /? in the quenched (a) and four flavours model (b). 

Fig. 4. v.e.v. of the inverse operator given by eq. (10) against (3 in the 

quenched case, 8"^, 10^ and 12"^ lattices. 
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Table caption 

Table I. Critical couplings and energies for different lattice sizes and flavour 
numbers. The superscripts x ^'^^ S denote values extracted from the chiral 
susceptibility and effective fermion action respectively. 
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lattice 


10^ 


10^ 


10^ 


8^ 


8^ 


8^ 


# flav. 


4 


2 





4 


2 





EX 


1.028(8) 


1.027(7) 


0.977(12) 


1.034(11) 


1.038(10) 


0.974(8) 


(3} 


0.2026(21) 


0.2231(18) 


0.256(3) 


0.202(3) 


0.221(2) 


0.257(2) 


E! 


1.02(2) 


1.030(15) 


1.030(10) 








/3f 


0.205(5) 


0.222(3) 


0.243(2) 









Table I 



11 



